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Abstract 

A system of generalized kinetic equations for the distribution functions of two-dimensional Dirac 
fermions scattered by impurities is derived in the Born approximation with respect to short-range 
impurity potential. It is proven that the conductivity following from classical Boltzmann equation 
picture, where electrons or holes have scattering amplitude reduced due chirality, is justified except 
for an exponentially narrow range of chemical potential near the conical point. When in this range, 
creation of infinite number of electron-hole pairs related to quasi-relativistic nature of electrons in 
graphene results in a renormalization of minimal conductivity as compared to the Boltzmann term 
and logarithmic corrections in the conductivity similar to the Kondo effect. 
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Introduction 



Recent discovery of two-dimensional (2D) allotrope of carbon, graphene, and experimen- 
tal demonstration of its massless Dirac energy spectrum has initiated a huge experimental 
and theoretical activity in the field (for review, see Refs One of the most interest- 

ing aspects of the graphene physics from theoretical point of view is a deep and fruitful 
relation with the quantum electrodynamics.-'^'^'^i^'^i^ In particular, anomalous transport 
properties of 2D Dirac fermions, such as finite conductivity of order of e'^/h in the limit 
of zero charge carrier concentrationii'iSi'i^iii^ii^ii^ can be associated with a specific quantum 
relativistic phenomenon known as Zitterbewegung.- The current operator of non-relativistic 
electron commutes with its kinetic-energy Hamiltonian and does not commute with the 
potential-energy one. Yet, it is vice versa for the Dirac electrons that is a reason for the 
Zitterbewegung. The same commutation properties hold for graphene in the case where the 
potential does not cause Umklapp process. Qualitatively, an impurity potential acting on 
non-relativistic electron creates random friction-like force which causes finite conductivity. 
This is expressed quantitatively in the standard theory of electronic transport in disor- 
dered metals and semiconductorsi^'i^'i^ii^S by deriving and solving the classical Boltzmann 
equation. The impurity potential action on the Dirac electron can not be described within 
such a simple picture. Despite this important difference many authors exploited the classi- 
cal Boltzmann equation to analyze electron transport in graphene.— ii^i^'2^i22i>2i Rigorously 
speaking, it is not clear what will be the limits of its applicability in this unusual situa- 
tion. Our work presents a consequent derivation of kinetic equations for the 2D massless 
Dirac fermions. Some of our results for the static conductivity are similar to those obtained 
by various quantum-field theory methods.— 'i^'i^'^'i^'i^ The approach based on the kinetic 
equation provides an alternative view on the anomalous transport properties of graphene. It 
can be easier generalized for more complicated situations such as strong electric fields, hot 
electrons, etc. These issues are beyond the scope of the present work. We will not consider 
also the effects of Anderson localization and antilocalization2^i2^i2^i^ in graphene restricting 
ourselves by the case of a weak disorder in the leading-order approximation. As we will 
see even in this case the problem turns out to be very nontrivial and instructive. We will 
prove that for not too small doping the standard Boltzmann equation with the scattering 
amplitude specific for massless fermions does give the leading term in the conductivity and 
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will find corrections to it due to the Zitterbewegung. In particular, these corrections have 
an interesting temperature dependence similar to the Kondo effect. 

A general idea of the approach used here is traced back to seminal papers by Kohn and 
Luttinger.^^ Starting from Schrodinger equation for noninteracting electrons in a random 
impurity potential they consequently derived the kinetic equation for diagonal (in momen- 
tum representation) matrix elements of the one-electron density matrix in the cases of weak 
potential or small impurity concentration. In these cases the kinetic equation turned out to 
be identical with the classical Boltzmann equation. Even for the simplest system to which 
it was initially applied, the Kohn and Luttinger treatment^^ proved not simple. For multi- 
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component systems one may also follow the route of Ref|! 
system for distribution functions in the momentum space - usual ones and functions that 
describe inter-subsystems transitions - but complexity of deriving such kinetic equations 
sharply increases. 

Several established formalisms exist nowadays, which automate the above derivation as- 
suming existence of some kinetic equations in principle. A partial list includes Kadanoff- 
Bayn^, Keldysb^i^, Zubarev nonequilibrium statistical operator (NSO)^ (for the NSO 
method, see also recent reviews^*^) and Peletminskii-Yatsenkc^^ methods. The Keldysh, 
NSO and Peletminskii-Yatsenko methods have close rationales. Namely, existence of an 
asymptotic density matrix which allows for Wick-rule decoupling of the creation and an- 
nihilation operators product averages is assumed in these methods. The consideration of 
non-equilibrium at strong interactions benefits using the Keldysh method which is distin- 
guished for highly developed diagram technique. At weak interactions, however, when the 
Born approximation is applicable the simplest approach in our opinion is with the NSO 
and Peletminskii-Yatsenko methods. This is because in the Born approximation, closed 
equations for the averages of gross variables, generalized kinetic equations (GKE), which 
describe non-equilibrium of interest (provided that such variables are declared in advance) 
were derived within these frameworks in late 60's once for all.— "^"^"^ 

In this paper we obtain and asymptotically solve GKE for spatially homogeneous 
graphene in order to calculate the linear-response conductivity. The main difference with 
the canonical case^^ is that for graphene the diagonal in the momentum representation av- 
erage density matrix is still two by two matrix in the pseudospin space, its off-diagonal 
elements describing the Zitterbewegung. This makes the GKE structure, on the whole, es- 
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sentially different from that of the classical kinetic equation. The structure of the paper 
is the following. In section I we present original expressions for the Hamiltonian, current 
and coordinate operators. In section II we specify the gross variables appropriate for the 
kinetics in spatially homogeneous case, which in fact are all the density matrix elements, 
and on their base concretize GKE regarding the interaction with arbitrary static impurities 
as a perturbation. In section III, assuming presence of a thermostat, we consider linear 
response regime (the case of small electric field) and express the linear static conductivity 
via two unknown functions of the one-electron energy. These functions, together with a 
subsidiary function of energy, satisfy a coupled system of linear integral equations resulting 
from linearizing GKE in electric field strength. In section IV, we solve the linearized system 
and calculate the conductivity within an ultraviolet cut off Dirac-delta impurity potential 
asymptotically in a controllable small parameter, using methods of solving singular integral 
equations. In section V we discuss the results obtained. 



I. PRELIMINARIES 

We proceed with the Hamiltonian of two-dimensional massless Dirac fermions describing 
charge carriers in graphene if one neglects the Umklapp processes between valleys K and K' 

n, = vY,%{'r.p)m^ (1) 

p 

where p is the momentum vector, v is the velocity, 

= ( <i <2 ) (2) 



^pl 

V^p2 

are two-component pseudospinor operators, 1,2 labelling the sublattices, and 



10/ I i 



(3) 



are the Pauli matrices in the pseudospin space. We will neglect here real spin and valley 
indices. The Hamiltonian ([1]) can be diagonalized using the unitary transformation matrix— 



where (f)p is the polar angle of the vector p. Hence the new electron operators given by 

1 / ^pi + e-^'^p^'pa 



V 2 \^ Vpi - e-*'^pVp2 
and 

-i>=^i^p=(4i ^u) (6) 

are the annihilation and creation operators of the conduction and valence band electrons. 
Thus we have 

*P = ^pSp, K = ^lUl (7) 

and 

Ho^v J2 ^P^s^ps- (8) 

p,s=±l 

In what follows we will consider the simplest case where electrons experience action of a 
scalar potential V (r) presents. The interaction Hamiltonian in this case is given by 

^int ^S-'J2^{p- p') *t,*p, = J2 ^IVPP'^P' (9) 
pp' pp' 

where 5" is the graphene layer surface area, V (q) is the Fourier transform of V (r) and 

Vp„' = -V {p - p') .a , J ■ (10) 

The current density operator in the new variables reads 

J = e^;5]*t^^p = J]stjpSp, (11) 
p p 

where 

jp = evUlrUp = {jpxjpy) , (12) 

and 

(cos (f)p —i sin 0p \ / sin 0p i cos 0p \ 

, Jpy = . (13) 

i sin (f)p — cos (t>p j ~ / 
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Off-diagonal elements of the current operator correspond to the Zitterbewegung processes.- 
For the x and y-components of the current, we further obtain 



cos 0r. —I sm I 



I sm 0p — cos 0p 



^ [cos0p (4,i^p.i - 4,-i^p.-i) - ^sin0p (^4,i^p.-i ~ 4,-i^p.i 



(14) 



and 



sm 0p z cos 0p \ _ 



p 

XI P^'^P (4,l^P:l -4,-l^P.-l) +«COS0p (4,1^P.-1 - 4,-l^P.l 



-z COS </)p — sm (/)p 



ev 



(15) 



At last, the electron coordinate operator which is necessary to derive the field term in 
the kinetic equation reads 

R = ^J]M/t.Vxl/p, (16) 



where V is the gradient operator with respect to the momentum p. Using the above unitary 
transformation, we get 



R = ^ E -.K^ (^p-p) = ^ E -l^-p + ' E -.K 
p p p 

= E \' (4i^^pi + 4,-i^^p.-i) - \ (4i^pi + 4,-i^p,-i 
p 

+5E'^'*p(«p.«p.-i+?i,-i«pi)- 



(17) 



To simplify R let us perform additional gauge transformation ^p^ — > e~*2<^p^p^ that re- 
tains the Hamiltonian Tio unchanged but renormalizes the coordinate and current density 
operators to take the form 



^ E (^4mV^P- + ^Cmep,-mV0p j = Rintra + Ri, 

p,m=±l ^ ^ 



and 



J ^ ev m ( — ^pm'Cpm + iipmip-mP^ 4>p j — Jintra + Jinter, 

p,m=±l ^ 



(18) 



(19) 
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respectively. Here we have separated explicitly intraband (electron-electron and hole-hole) 
and interband (electron-hole) contributions. Note that 



4eR = J, 
dt 

as it should be. The interaction matrix elements are thus transformed to 

v'pp'^v^(p-p') 



(20) 



COS I sm ^ P 



I sm — ^-E- cos 



(21) 



II. THE BORN-APPROXIMATION KINETIC EQUATIONS 



A. General outline 



The basic idea of the methods of Refs ]33lJ36l is a concept of so called "coarse-grained" 
dynamics. To apply the formalism we are, as noted in Introduction, to suggest the gross 
variables P, averages of which (P) at the kinetic stage of the evolution are believed to satisfy 
GKE. It was proven by Kohn and Luttinger— that, if V (r) is due to random impurities, the 
diagonal elements of the one-electron density matrix in the momentum representation aver- 
aged over weakly perturbed non-equilibrium ensemble are self averaging over the impurity 
configurations and do obey such a reduced description, at least for weak enough potential 
or small impurity concentration. Our problem is formally different from standard one only 
in existence of the interband operators. Therefore we choose the following gross variables 



^p-i^P-i 
^pi^p -1 

V 4,-iepi I 



(22) 



the components of this vector being the second-quantization form of the above matrix 
elements. The corresponding "quasi-equilibrium" or "coarse-grained" statistical operator 
(QSO)^^ is given by 



Pq = e 



-'f-E(^pi4i5pi+^p24.-i«p.-i+^p34i«p.-i+^p*34.-i«pO 



where 



$ = In Tr 



■ - E(^pi4i«pi+^p24,-i«p.-i+^p34i«P'-i+^p34,-i«pO 

e p 



(23) 



(24) 
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is the generalized Masseu-Plank function. As at the equihbrium, are parameters conju- 
gated to Pp in the sense that^'^'^'^ 

(5$ 



(25) 



This QSO is second-quantization representation form of the general density matrix in the 
which allows for Wick rules. Following from it explicit connection between (Pp) and Fp, 
however, bears no new information. 

To obtain GKE one uses the closure condition {Pp)^ = {Pp) assumed^i^ only for the 
gross variables, which results in 



^^{p)^ = ^{[no + n■,^uP]) 



(26) 



where the averaging in the right-hand sides is performed over NSO obtained from QSO via 
an explicit formal prescription—!^. In our case this averaging is also to incorporate one 
over the impurity configurations. Note that in all known cases with weak interaction the 
operators Pp obey closed microscopic dynamics with the unperturbed Hamiltonian 



(27) 



where cUpq is a known matrix. It can be shown that in our case Eq. (1271) sustains even if Tio 
includes the interaction with an electric field E along the x-axis 



n 



ef 



-eE.R=-eE ( ^.V^p™ - 

p,ni=±l 



p,-m 



(28) 



At that occurrence, the matrix tUpq contains linear in E off-diagonal elements, some of which 
involve the gradients of the momentum-conservation delta function. 

In the second-order approximation with respect to Tii^t Eq- (126|) can be transformed^S'^'^ 
to GKE, which have the following form common to all applications 

^=.5^.pq(Pq) + J« + j(^), (29) 
q 

where the generalized collision integrals of the first and second orders are given by 



and 



Jr 



(2) 



lim 



5 7^ 

n,^,(t),[Pp,n,r.,] + tJ2f^Pi 



(1) 



(30) 



(31) 



respectively. Here the dependence of the interaction Hamiltonian on time t is according 
to Heisenberg picture with Tio which, in addition to the kinetic energy, may include Tief- 
The term in Eq. fl3ip . which involves Jp \ leads to cancelling possible contributions that 
diverge in the thermodynamic limit S ^ oo out of Jp \ This is fair analog (in the Born 
approximation) to "connected-diagrams" statement in the diagram techniques.— 1^>^ 



B. Average current density 

Using Eq. f[T5]) the average current density can be expressed via the basic averages (Pp) 
as follows 

lj\ ey r 

jx = = J [(/pi - /p,-i) COS 0p + 2 Im (^pi) sin 0p] (fp, (32) 

where, by definition, 

fps = (^ps^ps) ) 9ps = {^ps^P,-s) = fi'p.-s- (33) 

Let us now introduce the electron and hole distribution functions 

np = fpi, Pp = l- fp,~u (34) 

which are, of course, real, and "anomalous" distribution function Qp = Qpi, which is complex 
in general. In the terms of these functions Eq. ( l32l) is written as follows 

ev 



i2nf 

ev 

(2^ 



[(^q + Pq) COS 0q + 2 Im (g^) sin 0q] d'^q 
j [A^q cos </)q + 2 Im (5(q) sin 0q] d^g. (35) 



C. Derivation details 



Let us now specify Eqs. 0291) - (1311) for graphene. To this end, consider all the prerequisites 
of the calculations required. Using Eq. fl^ . we find 



COS I sm ' 



1 ^ — ^ , + / >jWO 7. b Dill 7^ 1 

W,. = S-'5:V(1-1')S; , ^ ^^^^^ - E„ 



11' 



% sm „ ' cos ' 



^ 2^ r (1 - I ) ( cos ^ ^x^M',m + sm ^ ^iV^',-™ 1 

11' ™ LI ^ / 



(36) 



ir,m=±i 
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and so 



^0 

J — oo 



ir,m=±i 



cos 



pi-<Pi' 



e + im [ti - tv 



I sm 



e + im (e; + e;/) 



1',— m 



(37) 



Consider further the commutators of the gross- variable operators with Tio- By a straight- 
forward calculation we obtain 



I,n=±l ^ 



djcju eE_ 



psSps 



2p 



y sm </)p 



and 



[^0,^p/p-s] 



X ( nt;p^J„^q„ - ieE^l^^^ + ^ sin 0q4„^q,-n ) , ^i^^p. 

.q,n=±l ^ / 

Hence the "precession" terms in the right-hand side of GKE, see Eq. fl29|) . are 
^2^^pw(ei,/qs) = -e^TT- + Im^p,, 



sm I 



q 



ps,q— s 



{(Isi^s) = 2isvpgps - eE 



dgps 
dpx 



, eE sin i 



ifps fp,~s) 



(38) 



where we have used the notations introduced in Eq. (!33ll . 

To calculate collision integrals for the "normal" distribution functions fps and the anoma- 
lous ones Qps we are to perform the commutation twice - first time to commute the gross- 
variables operators with Hint to obtain J7p^'' and the second to commute the result of the 

(2) 

first commutation with £, to obtain J'p . Following this route we get 



+ i sin 



qq 

0q 



COS 



^P,q "P,qV SqsSq's 



^qs^q 



^P:q^qs^q',-s ^p,q'^q,_s^q' 



(39) 



2 \^-F:q^qsSq 

Averaging this expression over QSO gives J7p^'' [fs] = 0. Then, performing the second com- 
mutation using Eqs. (!37|) and (!39l) . after straightforward calculations we obtain the second- 
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order collision integral for /, 



ps 



^>p - 0q) Re (^q) 



:7-f [fs] = S-'J2\V{p-ci)\'f^ssm 

— ifps - Us) + s sin (0p - 0q 



1 1 

+ 



^E) ~l~ ^0 '^I 



2 cos 



2 V^P 



(40) 



Further we have 



+ 2 sin 



qq 

0q 



COS ■ 



— (5pq - ^pq')^i.^q',-s 



2 ^^'-'pqSqsSq's <-'pq'Sq,-sSq',-s 



(41) 



Making the second commutation with the use of Eqs. fl37j) and fl4T|) . after straightforward 
calculations we obtain the second-order collision integral for Qps 



J^'^ [9s] = -S-'J2\V{p-ci)\'\2 

qq' 



COS 



2 V^P 



{9ps - Qas) 7r(5 (cp - 



9ps ~l" fi'qs 



- sin (0p - 0q) X 



(/qs ~l" /q,- 



7r5 (ep - 



+ 



s . 

- sm I 
2 



/qs - /q -s _ 2^_g g-j^2 0p " 0q 9ps + fl'qs 



Cp + Cg 



(42) 



2 Cp + Cg 

This equation describes the Zitterbewegung effects, that is, creation of electron-hole pairs 
during the charge carrier propagation. Note that in both Eq. (HOj) and (142!) the configura- 
tional average of the potential Fourier transforms squared is implied. 

Putting together Eqs.( |38ll . (H0ll and ( H2l) we can write the final set of GKE for the "normal" 
and "anomalous" distribution functions. It is more convenient, however, to transform these 
GKE to a system of equations for the functions 



Dp= ^ /ps - 1 = rip - pp, Np= ^ sfps + 1 = np 



(43) 



s=±l 



s=±l 



and for Qp. Let us remind that the latter and A^p define the average current, see Eq. fl32l) . 
We have for Dr, 



dt dp^ 



2tt 



f cos' ^^-^ {Dp - D^) 6 {e, - e,) , (44) 
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which does not involve gp at all, while the equation for A^p reads 



dNp dNp 2eEsin0p 

+ eE ^ Im Qp 



dt dpx 
_ 2ti 



P 



5^ 1^ (P - q)l' sin (0p - 0q) Re (g^) 



cos 



q 

2 0P 



- (A^p - ATq) + sin {<Pp - 0q) Im {g^ 



~l~ Cp Cq Cp 



^ (ep - eg' 



Finally, the equation for the complex function gp proves the following 



+ 2 cos 



q 

2 0P 



|2 ) « . 

--sm, 



M^q 



(^P - c/q) (ep - eg) + 



S (ep - eg) + - 



i 1 



TT Cp - tq 



+ 



1 iVq 

27r Cp + 



. . , , . 2tgp + g ^ ^(f)p- 

sm (0p - 0q) sm - 

TT ep + eq 2 



(45) 



(46) 



III. THE LINEAR RESPONSE REGIME IN ELECTRIC FIELD 



In general, Eqs. flH|) - fH^ are quite complicated. Further we will consider only the regime 
of linear response, that is, the case of weak electric field. We will restrict ourselves also by 
the case of stationary field and neglect its effect on the collision integral. It can be shown 

(2) 

that linear in E corrections to jTp restore some second order terms of the perturbation 
expansion of the exact field term considered in the Luttinger-Kohn formalisnt^, which lies 
out of our scope. 

At E = 0, like the classical kinetic equation, Eqs. (jUj) - (1461) have equilibrium solution. 
In our case it is three arbitrary function of the energy ep = vp. To develop meaningful 
linearization of GKE and the current for small E, we assume, following Ref. [291, presence of 
a thermostat which role is only to establish true equilibrium with a temperature T. Under 
this condition the equilibrium distribution functions nop and pop become, of course, the 
electrons and holes Fermi functions, respectively, with unique chemical potential /x. Now let 
us consider the linearization of GKE derived above in detail. 
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A. Linearized equation for Dp 



Replacing in the field term of Eq. (H^ the distribution functions by their equilibrium value 
and linearizing the corresponding collision integral, we obtain the following equation 

eE^^ = -% j:\Vip- cos^ iSDp - 6 (e, - e,) , (47) 

where 

^o(e) = ^E^-^±^' (48) 
Eq. (H7I) is quite similar to the classical kinetic equation and so is routinely solved exactly. 
The solution reads 

6Dp = -eEvT (e,) cos 0p, (49) 

T (ep) being standard elastic transport relaxation time given by 



^ |^$^|l^(p-q)|%in^(0p-0,)5(e 



p ^q) 



-N- 



(ZTTf ) 



2t , 



sin^ (t)d(t), (50) 



where A^imp = c/fi is the impurity concentration per the area unit, c and Vt being atomic 
fraction of impurities and the graphene crystal cell area, respectively, and ?7(|q|) is the 
Fourier transform of one-impurity potential. In deriving Eq. (1501) and what follows we 
adopted that c ^ 1. The factor sin^ (p instead of standard^^ one 1 — cos is a consequence of 
the chiral character of charge carriers which leads to the suppressions of back scatteringi^ii^. 



B. Linearized equations for A'^p and 



Replacing in the field terms of Eqs. (H5|) and (l46ll the distribution functions by their equi- 
librium value and linearizing the corresponding collision integrals, we arrive at the following 
system of coupled equations 



eE 



dNo (ep) 2eE sin 0q 
Im ^0 (cp) 



dpx 



q 



1 1 1 

+ 



ea + e„ IT ea - er 



Re6g^ 



cos 



,2 S^P 



{6Np - 6N^) + sin (0p - 0q) Im 6g^ 



S (cp - eg) 



(51) 
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and 



— 2iepOg-p + eE — - — — + i 



eE sin 1 



dpx ' " 2p 
2 I i 



2 ^ 



^ [iVo (ep) - 1] = 
i 1 



+ 2 



- 5^q) S (ep - eg) + 



S (ep - eg) + 
i 6gp + 6gq 



sin (0p - 0q) 



cos 



2 S^P 



+ 



1 SN^ 

271 €„ + e. 



sm 



TT eg - ep 
2i5gp + Sg*^ . 



TT ep + eg 



sm 



2 V^P 



(52) 



where 



Nn(e) 



+ 1 + 1 



(53) 



and g^ (ep) is an equilibrium "anomalous" distribution function satisfying the equation 



2«ep^o (ep) = 5^ 1^ (P - q) 



-0^0 (eg) 



^ (e« - eJ + 



TT ep - eg 



sm I 



+ 



2i 

TT 



9o lepj + go [tg) ^ (pp_ 

COS 



^ _ 9o (ep) + go (eg) ^.^2 0p " 0q 
Crj -|- e^ 2 



+ 



1 No (eg) 
2vr ep + eg 



sm 



Because 



/27T /'2n 
I^(P - q)l^sin(0p - 0q)(i0q oc J U (^a/p^ - 2pqcosx + q 



sin xdx = 0, 



the terms containing Dq (eg) and A^o (eg) in the above equation for go (ep) give zero contribu- 
tions. As a result, if even non-zero go (e) is purely real. In contrast with A^o (e) and Do (e), 
this real function has no influence on the non-equilibrium 6Np and 6gp and so drops out of 
our treatment. 

Let us search the functions SNp and 6gp in the form 



SNp = —eEv (ep) cos0p 



and 



(54) 



(55) 



5gp = eE'f (ep) sin</)p, 

respectively. Substituting these forms into expression for the average current (1351) yields the 
following expression for the conductivity 

„2 roo 



a 



2ttv 



Im7 (e) 



ede. 



(56) 
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Substitute further Eqs. (l5^ and (155!) into Eq. (l5T]) . Taking into account that Im^fQ (ep) = 0, 
we manage to show that all non-zeroing terms of the resulting equation contain common 
factor cos0p. Dividing by this factor we get the first linear-response GKE 



der. 



Tier 



oo 



KQ + P Q-P 



1 1 

+ 



Re'j{eg)dq, (57) 



where the kernel in the integral term is given by 



To (p, q) 



2-K 



U [\/ p"^ — 2pq cos + j 



sm 



It is expedient to note that 



To {p,p) 



1 



(58) 



(59) 



vrr (ep) 

Performing quite similar transformations in the "anomalous" kinetic equation, Eq. fl52p . we 
find that all non-zeroing terms of the resulting equation are proportional to sin0p. Can- 
celling this common factor, we obtain the second linear-response GKE 



2iep7 (cp 



1 - No (ep) dDo (ep) 



der. 



imp 



|t/(p-q) 



sm 



dDo (eg) VT {eg 



{2nY J '~ { 2 

+ [1 + cos (0q - 0p)] [7 (ep) (1 - COS (0q - 0p)) TtS - e^ 

7 (ep) + 7 (eg) cos (0q - 0p 



+ 



deg ep-tg tp + eg 



+ I 



eg - ep 



• n ^ M 7 (ep) + 7* (eg) cos (0q - 0p) , . 
-2 [1 - cos (0q - 0p)J — \ d q. 



er, + 



(60) 



C. Linear- response GKE in the energy variable 



In what follows we will make overall use of the energy variable e = vp. Let us summarize 
linear-response GK, i.e Eqs. ( 1571) and (1601) . in terms of e. We have 



de r(e) r (e) Jo \uj + e u - e 



^ + —^\Re-i(uj)duj (61) 



and 



- 2ie7 (e) 



IV 



1 - No (e) 9Do (e) 



7(e) 



re 



+ iA (e) 7 (e) + « 



9e 

r+ (e,c^ 
e — oj 



To {e,uj) 



r_(e, cj) 
-7 {uj) + ^^7 ^ 



V [uS) ^ dDo iui) VT (a;) 
e + oj doj e — ijj 

duj. 



duj 
(62) 
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where the kernels are given by 



Tn (e,uj) 



UJ 



{27rv) JO 



(27rf) JO 



27r 



2-K 



f/ ( — A/e^ — leuj cos </> + cj^ 



sm 



e Co" 



1 



U \ —-J — 2tuj cos + 0^2 

V 



[1 ± COS (p) COS (pdcf) (63) 



and, using these kernels, two-particle energy shift by 
Afe) = 



Fq (e,u;) +r+ (e,u;) _^ ro_(e, c^) - r_ (e, uj) 



duj. 



(64) 



e — UJ UJ + e 

Let us introduce dimensionless kernels (e, cj), a = 0, ± by the following identity rela- 
tions 

To (e, u;) = u;$o (e, u;) , T±{t,uj) = ±uj^±{t,uj) , (65) 

with which and Eq. fl63l) the formulas of these new kernels being straightforward. Note that, 
by construction of the kinetic equations, |$q (e, uj)\ <^1. Making use of the kernels $q, (e, uj) 
and of new functions of the energy defined by 

1 dNo (e) 



2/o (e) + 



7r$o (e, e) de 



e^{e)=vfie), / = /i + z/a 



(66) 



we arrive at a modified conductivity expression and two coupled singular integral equations 
for /o (e) and / (e). We have 



g2 roo 

a = aB + — [/2 (e) - /o (e)] c?e, 
2vr Jo 



where 



SATo (el 



{27Ty Jo '^'o(e,e) 9e 



(67) 



(68) 



is the classical Boltzmann conductivity^*^, with the Born impurity scattering cross-section 
modified due to chirality as noted above, in which of electrons and holes contribute additively, 
while the additional non-classical term is due to the Zitterbewegung. The integral equations 
are as follows. The first one is real and homogenous 

1 



n% (e, e) [f^ (e) + /o (e)] - H % (e, uj) 

Jo \^ + 



+ 



+ e uj — e 

and the second equation is complex and inhomogeneous 



/i (u;) duj = 



(69) 



[A (e) + ^7r% {e, e)] / (e) 



UJ — e 



UJ + e 



duj 



e + UJ 



fo{uj)duj = F{e), 



(70) 
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in which 



and F (e) = Fi (e) + iF2 (e) , where 



A (e) = 2 + e-^A (e 



l-iVo(g) I IdDo (6) 



2e 

$0 (e,a;) 



2 9e ' 



1 9A/'oH 1 9DoH 



uj + e did uj — e duj 



doj. 



Eq. flTOl) is, in turn, equivalent to two real equations for /i and /2 

(e, cj) $_ (e, 



A (e) A (e) - 7r$o (e, e) h (e) - /' 

^0 



UJ — e 



and 



7r$o (e, e) /i (e) +A (e) h (e) - 

$0 (e,o;) 



$+(e,a;) $-(e,u;) 



u; + e 



uo — t UJ + e 

fo (uj) duj = Fa (e) , 



/2 (uj) duj 



(71) 



(72) 



/i(a;)rfcu = Fi(e) (73) 



(74) 



respectively. 



IV. A SHORT-RANGE IMPURITY POTENTIAL MODEL 

In this section we solve the system of the integral equations derived above, see Eqs. (!69|) 
and flTOj) . or the real equivalents of the latter - Eqs. fl73l) and (1731) . and then calculate the 
conductivity using Eq. (1671) for a weak extremely short-range impurity potential. 



A. Formulation of the model 



Let us consider a zero-range impurity potential, which we define by U {r) = UoQ5{r), 
where the parameter Uq has dimension of energy. For this potential, U (q) = QUq so the 
kernels become independent of their energy arguments and all equal 

Here where we have restored the Planck constant h to stress that $ is dimensionless. This 
approximation poses no formal problem as regards the integral terms in Eqs. (!69|) and 
rOjl but A (e) acquires an ultraviolet logarithmic divergence. Thus some ultraviolet cut off 
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procedure should be introduced. We define (e, cu) = $ at < e < ec, < c<j < ec and zero 
otherwise. Then the simple calculation using Eqs. flTTj) and fj64|) yields 



A fe) = 2 - 4$ 



:duj 



$ln 



(76) 



at e < ec and A (e) =2 otherwise. This function has infinity breakpoint at e = ec, which 
has no physical meaning. Analysis of a general case with a finite-range potential shows that 
the model is reasonable assuming that we are interested only in small enough energies in 
comparison with ec = hv/ro where tq is a characteristic radius of the potential. On the other 
hand, the noted mathematical property of the model A (e) allows us to solve a part of the 
obtained singular integral equation exactly, see Appendix. 



B. The model integral equations and their solution 

Using the introduced model we obtain from the following system of singular integral 
equations for e < ec 

/2 (e) + /o (e) - - = 0' (77) 

A (e) /i (e) - n^f2 (e) - 2$ 4"^^^^ = (^) ' (78) 

Jo ^ — ^ 

TT^A (e) + A (e) h (e) - 2$e 4^^^ " ^ f = F, (e) , (79) 

Jq LO — e Jq uj + e 

where 

At e > ec we have 

A(e) = ^Fi(e)^-(4e)-\ (e) = 0, (81) 

so /2 (e) = in this range. Yet, no definite a priori information on /o (e) at e > ec can be 
deduced in the considered model. 

This set of singular integral equations can be solved using the well-developed methods'^- 
of complex calculus, which is presented in Appendix. Using Eq. (11181) from Appendix, we 



1 dNoiu) 1 dDo{uj) 



u + e duj 



id 



duo. 



(80) 
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get for the function which directly determines the Zitterbewegung conductivity in Eq. ([6 



/2 (e) - /o (e) 



2(7r<l>)^ 



(27r<I))' + A2 (e) vr$ 



+ 



A(e) 



A2 (e) 



X 



TT 



-e 



' Fi (u) udoj 



A2 (cu) 
2$A ie) 



(271$) 



2 tU^ 



A2 



+ (27r$)^ 

2 (vr^)^ /o (6) 
(27r<l>)' + A2 (e) 

'/o(u;) wciw 



x2 a;^ — e 



2 ■ 



(82) 



f(27r<l>)'-^o y'A2 (cj) + (27r$)^ 

In turn, /q (e) satisfies Eq. (11191) which is closed Fredholm like integral equation with a 
kernel non-singular at e = cj. As we will see immediately below, the first term in Eq. (!82l) is 
compatible to a Zitterbewegung one obtained from qualitative analysis of the Kubo formula 
for ideal Dirac fermions.- The last two terms give rise a novel Zitterbewegung contribution 
to 0", which results in post-leading corrections, O (1) at most, for $ — > (see Appendix) 
and so neglected here. Note that, though local in the energy, the first term in Eq. fl82l) is 
just a result of solving the singular integral equations system, Eqs. ( 1771) - (!79l) . and not of 
plain approach when /o (e) and all the singular integral terms of the system are neglected 
in advance. 

Using the adopted approximation in Eq.f l56l) yields for the Boltzmann conductivity part 
in the units /h the foUowing^^ 



1 



0"B 



dNo (e) 
de 



de 



A^o(O) 

277$ 



27r$' 



(83) 



up to the terms 0(Tmin(/i^^, e^^). Let us now integrate the first term in right-hand side of 
Eq. fl82p over e using integration by parts. Thus we obtain the Zitterbewegung contribution 
to 0" in a pseudo-Boltzmann form. Using Eq. ( 172|) . we have 



(7r<I>)" 

O'ZB = O'B / o 

/o (27r<I>)VA2(e) 



dDo (e) 



, ^ dNo je) 
erzB (e) —^^ e 



e de 
dtTzB (e) dDo (e 



de 



where, by the definition. 



TZB (e) 



Ob 
e 



de 



de 



de. 



(84) 



du 



{27T^y + A2 (io) UJ 

is an effective Zitterbewegung relaxation time. 



(85) 
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C. Applicability of classical Boltzmann equation and analysis of the /i = case 



From Eq. (176|) we find at e <^ 



e 



A (e) ^ 2 (^1 + 2$ In - J (86) 
i.e. A (e) zeroes at the energy 

ex = e,e-^ = e.e''^'^^ (87) 

which is a striking analog of the Kondo energy scale in a problem of magnetic impurity in 
metals.-^ Existence of this exponentially small energy scale in the problem under consider- 
ation was established first in Refs Jll|]l2l . Estimations of contributions from Eqs. fl8^ and 



Eqs. (l85l) as well the neglected Zitterbewegung terms show that at 

ec > > max (ex, T) (88) 

the corrections to the Bornian conductivity fl83l) are at most finite in the limit $ ^ and 
thus can be neglected in comparison with ctb- This justifies using the classical Boltzmann 
equation for graphene, except the case of extremally small doping. 

Formally speaking, application of the theory developed here to the case of zero doping 
is doubtful. For example, the self-consistent Born approximation^- gives for this case dras- 
tically different results in comparison with the Born approximation. At the same time, our 
approach is formally exact in a sense of perturbation theory at $ — » 0. We will see that, 
actually, the classical Born-approximation Boltzmann equation does not take into account 
properly all terms of order of a part of such results from Zitterbewegung. 

Let us now perform integration in Eq. fl5^ assuming validity of Eq. flHB]) . which is fairly 
justified at $ ^ 0. This yields for the case /i = 



erzB (e) = 



1 

arctan arctan 

$7r 



(89) 



Substituting Eq. ( l89l) into Eq. ( l84l) . we obtain the integral formula for the conductivity in 
the undoped graphene {fi = 0) 



1 1 ^ 

a = (Tb <l + - 



TT /2, 2T 

— arctan — In — x 



dx 



cosh^ X 



l2 V'^ 

Given $ <^ 1, this formula has the following asymptotic behavior with respect to T 



(90) 



l + ^ + g(lnf)-\ T«6K 
S . .2 /, r \ -1 _ (91) 



1 + 16 1^4) ,T»eK. 
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It is seen that the Zitterbewegung correction at T = has the same order of magnitude as the 
Bornian conductivity and numerically even larger than it by a factor 7r^/4. The temperature- 
dependent corrections are reminiscent to those in the early theories of the Kondo effect by 
Abrikosov and Hamann.'^^ Thus we obtain the following conductivity ratio 

o-(oo) 4 ^ ^ 

V. CONCLUSIONS 

We have derived the second-order perturbational GKE for 2D massless Dirac fermions 
in graphene scattered by scalar impurity potential. We considered the GKE solution in the 
Dirac-delta potential model with the ultraviolet energy cutoff. Our principal result is the 
criterion given by Eq. fl88|) . which justifies using the classical Boltzmann equation, except 
for exponentially narrow interval of chemical potential and temperature. Our approach 
clearly demonstrated that the problem of conductivity at zero doping it fairly similar to the 
Kondo problem. In this case, we obtained the temperature dependent conductivity formula, 
Eq. fp}]) . which interpolates well between high-temperature (T ^ ex) and low-temperature 
(T <^ ex) ranges and remains finite at T = ex. Thus consistent asymptotic solving the 
integral equations that result from the derived GKE, we performed in this paper, proves 
equivalent to a partial summation of the perturbation terms oc In T. Similar procedure was 
carried out for canonical Kondo model using the NSO method in Ref. |39. 



By the noted analogy with the Kondo problem,"2S Eq. ( l90l) at T > ex is asymptotically 
correct in the controllable small parameter while we may not pretend to describe by 
it the low-temperature properties, in particular a (0) in detail. Nevertheless for all T, our 
kinetic equations by construction are more general than the Boltzmann one (even with 
scattering rate modified due to chirality of the current carriers). Therefore, the discrepancy 
by factor ~ 3.5 in the values of a (0) obtained, see Eq. fl^^ . makes the Boltzmann equation 
probably not very good starting point for generalizations, such as the self-consistent Born 
approximation.— Would the Kondo analogy goes pretty far, the observedi cr (0) ~ e^/h 
might be an evidence that the system enters at T ^ ex a non-perturbational strong effective 
coupling regime, where the conductivity attains so called unitary limit,— rather than the 



result of strong bare coupling $ ^ 1 adopted in Ref. 



16. 
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Appendix 

Here we present the solution of integral equations (!77|) - (!79|) . Subtracting Eq. (177|) from 
Eq. flTSl) we obtain purely algebraic linear equation 



A (e) A (e) - 27r<l>/2 (e) - 7r$/o (e) = F, (e) 



(93) 



Let us exclude now /i (e) expressing it via /q (e) and /2 (e). To this aim, we use Eq. (177|) 
in the form of problem of inverting the Cauchy integral 

2 

TT Jo 



4^^du; = fo (e) + f, (e) = p, (e) . 



(94) 



In what follows we make use of the celebrated Poincare-Bertrand permutation formula:^ for 
the Cauchy-type integrals along a contour C in complex plane 



c 



T — t Jc UJ — T 



dr 



c 



c 



V; (e, r) y (r, g;) 
(ti -t){u- t) 



dr 



duj — Ti'^ip (e, e) if (e, e) . (95) 



Putting in Eq. ([91]) e r, = if (r, oj), multiplying it by 

V;(e,r) 1 



r — e 



(r2 - 62) v/if372 
and integrating over r with the use Eq. fl95p . we obtain 



"Co 



Pi (r) dr 



vr/i (e) 



+ - 



c / *=c 



1 



dr 



I (r2 - e2) v^if^ Sev^if^ ' vr Jo 

7r/i(e) 2 I {lo') - I {e') dr h{u 



r2 _ e2) - r^) ,/i^372 



/i (cj) curfcu 



2e^ 



.2 TT 



rUjduJ, 



where (the integral below is the principal-value one) 



dr 



dx 



dt 



- e2) A/e2 - t"^ cos^ x - Jo 



t2 
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from which we deduce that general solution is 



/. (e) ^ -'-^^^ r + (96) 

provided that pi (e) is supposed to be known and Ci is an arbitrary constant. Consider now 
the behavior of the above solution at the interval ends. We have identically 

h (e) 




Thus, if pi (0) and pi {e^) are finite, /i (e) at the ends diverges if Ci 7^ and is zero if Ci = 
since = 0. Choosing Ci = 0, we obtain from Eqs. fl93l) and fl96|) for the function 

„M = ^l|l= (97) 

the following singular integral equation 

2 udu ^ ^ , , 7r$/o (e) - Fi (e) 



A (e) - / qi {uj) — + 27r$gi (e) 

TT Jo a; — e 



= p2 (e) . (98) 

Introducing new variables and functions by 

Vx = e, y/y = uj, qi (x) = qi {^/x) , p2 (x) = p2 (Vx) , (99) 

we convert Eq. fl98l) to the standard singular integral equation^I 

A (Vx) - I ^^-^dy + 2Ti^qi (x) = (x) , < x < x^ = e^, (100) 
Jo y - a; 

assuming ^2 (x) is known. Following the standard procedure^, we define the function of 
complex variable 

Q. (^) = ^ r ^rf!/. (101) 

which is analytic in the plane with the cut along (0,Xc) and 

lim Qi {z) = 0. (102) 

|2|^0O 
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The relations at z ^ x ±0 



qi {y) 



dy = Qt (x) + Qi (x) 



-^i Jo y-x 

map our equation onto the Riemann-Hilbert boundary value problem 

lA (v^) [Qt (x) + gr (x)] + 2vr$ [Qt (x) - Q^ (x)] = p2 (x) 



(103) 



or 



where 



Qt {x)-Gi (x)gr (x 



P2 [X) 



Gi ix) 



27r$ + iK {^/x) ' 
27r<l> - iA (v/x) 



(104) 



(105) 



27r$ + iK (v/x) ■ 

Proceeding, we are to solve the homogeneous Riemann-Hilbert problem of searching a regular 
analytic function VLi [z) satisfying 



nt (x) = Gi (x)(]7 (x), 



(106) 



which is considered below. To obtain the solution of the inhomogeneous problem following 
Refj37l we divide Eq. fll04p by Vtt i.x) and using Eq. fll06p obtain 

Qt {x) _ Q^ (x) _ _l P2 {x) 

Q.t {x) VL-{ (x) ~ VLt (x) 271$ + iK (v^) ' 

from which it immediately follows that 



Qi {z) = {z) 



P2 (y) 



dy 



(107) 



2m Jo nt (y) 27r$ + lA {^) y - 

where Pi (z) is an analytic function in the whole plane except may be points z = and 
z = Xc- The values of Qi [z) on real axis allows one to obtain gi (x) and its Cauchy integral 
using first of Eq. fll03p along with Eq. (11070 as follows 

1 1 p2{y) dy 



gi (x) = [nt (x) - (x)] 



+- 



+ 



• nr_(x) 

Qtix) 

A (y/j) Qt jx) 
27r$ - zA (y/x) 

27r<l>p2 (x) 
(27r$)' + A2 ( Vi) 



27ri Jq Qt (y) 27r$ + iA {^) y-x 
P2 (x) 



+ Pi (x) 



27r$ + iA (y/x) 



TC 



1 



P2{y) 



dy 



Qt (y) 27r$ + iA (yy) y-x 



2iPi (x) 
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and 



Jo y-x 



P2 (y) dy 



+ 2 



fir ix) 



2n Jo nf (y) 27r$ + %A (^) y-x 
P2 (x) 



+ iPi (x) 



27r$fi+ (x) 
" 27v^ - iA (v^) 

A {^/x) % {x) 



27r$ + ik {^/x) 

1 1 



P2 {y) dy 



n Jo nt (y) 2tt^ + ik {^) y 



X 



+ 2iPi (x) 



(27r<l>)' + A2 (v^) ■ 

Returning to the relevant homogeneous problem we assume that Qi (z) 7^ 0, 00 at 2; 7^ 0, 
Xc- Thus we arrive at the inhomogeneous problem for Infii {z) 

A(v^) 



InQf (x) — InQi (x) = \nGi (x) = — 2zarctan- 



27r$ 



Note that the end-point conditions are 



lim InGi {x) = in, lim InGi (x) = —in, 

x^0+ x— >Xc— 



;io8) 



where the first limit holds in general case and the second is the model property. Consider 
the following Cauchy integral^^ 

J_ r^\nGi{x) 
2m Jo 



UAz) 



-dx 



1 A(v/i) dx 

-— / arctan . 

7r Jq 27r$ X — z 



X — z 

This function satisfies lim|^|^o Ui {z) = 0. It has two regular nodes at the points z = and 
z = Xc, which is shown using integration by parts 

Xc-O 



TT / \ 1 1 / \ A (V^) 

Ui [z] = m (X — z) arctan 



I r\ , , d A(v/i) , 

H — / m (x — 2;) — arctan ax 

TT Jo ^ ' dx 27r$ 



--\n\{-z){xc- z)]+2(^ I — ^ ^-^^dx. 



(109) 



Thus the function e^'^^^^ can be taken for ^li(z), which satisfies limi^i^oo (-2) = 1 and 
hence it should sustain limui^oo Pi {z) = 0. For this Qi {z) we have on real axis 



nf (x) = e 



Ui{x±iO) 



■ — fr\ Syrctctn — 1^ — T 

TT JO 27rq> y — x^xO ^ g 



'S(x)^i arctan ^^^^ 



27r$ =F iA {^/x) 
A2 {^) + (27r<l>)^ 



.61 (x) 



(110) 
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where 

e.(.) = -i^arctanM^|)^. (ni) 
^ ^ TiJo 271^ y-x ^ ^ 

This yields with Pi = 

A(v^)e®i(^) 1 r= e-®i(2/)^2 (y) 



A e«i^-^ 1 f 
. /a2 ^ /:;:^ , ^o^^^2 tt Jo 



(27r$)' + A2 (v^) 



+ : 2 , (112) 



and 

1 r^Qiiy), 27r$e®i(^) 1 r= e-®i(f)p2 (y) (iy 



1 r"Qiiy)^y_ 27r$e«iw i /-^ 



/A2(Vi) + (27r$)'''-^o y^A2 (^) + (27r$)'^ ^ 

A {\/x) P2 (x) 
(27r$)^ + A2 (^) 

Returning to the energy variables, we obtain 



Z . (113) 



^ 2(7r$)V3(6) _ 2$A(6) r= e®(-^)-®(-^)/3(c.) a;da; 

(27r$)VAM6) ^A2 (,) + (2^$)2 io (^) + (2^$)^ - 



and 



TT^A (6) /3 (6) TT (2$)^ ... ^e(e^)-e(.^) 

^^^'^ (2.$)^ + AMe) ^A2(,) + (2,$)2yo (,) + (2.$)^ - 



where 



and 



/3 (6) = /O (6) - (116) 



e(e2) =ei (e2)+lnVe' (e^ - a;^) = 2$ T ^^^^^^^^ . In |x - e^l 



A^^) + {2TT<l>y 
ln|a; — e^l dx 

(l + $ln^)' + (7r$)'^(''~^) 



In 

oo 



+ 1 e2 



(117) 



($-1 + «) +7r' 



the constant $ In being omited since only difference © (e ) — © {uj ) enters all formulas. 
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/^From Eq. (11141) and the definition of pi (e) we obtain tlie connection between two 
functions of interest 

\2 



/2(e) 



F^{e) 2(7r$)^ + A2(e^ 



7r<l> (27r<l>)^ + A2(e 
2<I>A (e) 



v2 tu^ — 



(118) 



(e) + (27r$)^ "^0 ^ (cu) + (27r<l>)' 

Using Eqs. (II 151) . (11 181) and (I116P along with Eq. (!79!) and the Poincare-Bertrand formula, see 
Eq.( l95l) . we obtain closed integral equation for the function /o (e). The equation reads 

1 



A(6)/o(e)-$ 



uj + e 



/o {uj) duj = F3 (e) 



where the kernel and inhomogeneity term are given by 



(119) 



UJ 



[{2Tr^f + 2A2 (e)] e 



e2 -e 



A2 (e) + (27r$)%/A2 (u) + (27r<l>)^ 



e [(27r$)^ + 2A2 {u 
(27r$)' + A2 (w) 



and 



(e) = -F2 (e) 



7r 



2e 



C(j2 — e2 



+ Q{e,uj) 



respectively, while 
Kie,cu) 



')A(r) dr 



y^A2 (uj) + (27r$)' (^2 _ ^2) (^2 _ ^2) y/A2 (^) + (2vr$)' 



(120) 



(121) 



(122) 



is another kernel. Note that, like K {e,uj), the kernel Q {e,uj) is non-singular at = e 7^ 0. 
Also we have 



lim K (e, uj) = 2euj 



dr 



[r2 - e2) (^2 _ ^2) 



e2 — Ci;2 



and 



lim Q (e, uo) = 2 



UJ — e 



' ' ■ ' co'2 — e2 e + cij 



Further, the leading terms in $ at $ ^ 1 of the inhomogeneity term is 



lim f3 (e) = -F, (e) + 1 /"• ^lM|:<i. = (.) . 



(123) 



(124) 



(125) 
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It is clearly seen from the above equations that /o (e) = O (1) at $ — > 0, which results in 
the conductivity formula obtained in the main text. 
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